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Abstract.  The  recrystallization  behavior  of  cold-rolled,  commercial-purity  titanium  was  studied 
experimentally  and  with  Monte-Carlo  (MC)  modeling.  Utilization  of  EBSD-OIM  as  input  for  MC 
modeling  resulted  in  realistic  predictions  of  recrystallization  kinetics,  microstructure  and  texture, 
which  were  in  good  agreement  with  experimental  results.  MC  modeling  of  recrystallization  kinetics 
predicted  that  non-uniform  stored  energy  distribution,  heterogeneous  nucleation  of  recrystallization 
and  recovery  in  combination  leads  to  a  negative  deviation  from  linear  JMAK  kinetics.  It  was  found 
that  concurrent  recovery  that  takes  place  during  recrystallization  is  an  important  process  that  controls 
both  the  overall  recrystallization  kinetics  and  the  deviation  of  linear  JMAK  kinetics.  On  the  other 
hand,  non-uniformly  distributed  stored  energy  itself  has  little  effect  on  the  negative  deviation  from 
JMAK  kinetics  but  intensifies  the  deviation  when  heterogeneous  nucleation  is  combined.  Modeling 
results  also  revealed  that  heterogeneous  nucleation  of  recrystallized  grains  and  their  early 
impingement  in  local  areas  of  high  deformation  are  essential  for  producing  a  log-normal  distribution 
of  grain  size  and  a  typical  recrystallization  texture  in  rolled  titanium. 

Introduction 

Recrystallization  is  of  considerable  theoretical  interest  and  practical  importance  and  thus  various 
approaches,  such  as  analytical,  cellular-automaton,  Avrami,  and  Monte-Carlo  (MC)  techniques  [1-4], 
have  been  developed  and  utilized  to  describe  recrystallization  phenomena.  Among  those  methods,  the 
MC  approach  is  very  attractive  because  it  can  describe  a  number  of  different  aspects  of 
recrystallization,  including  the  overall  kinetics  and  the  evolution  of  microstructure  and  texture. 
However,  detailed  quantitative  information,  including  initial  grain  morphology,  crystal  orientation, 
and  the  distribution  of  stored  energy,  is  needed  for  MC  modeling.  Unfortunately,  unrealistic 
assumptions,  for  example,  a  random  distribution  of  the  nuclei,  homogeneous  distribution  of  stored 
energy,  or  random  starting  texture,  have  been  employed.  Therefore,  predictions  of  recrystallization 
phenomena  based  on  such  unrealistic  assumptions  may  lead  to  marked  deviations  from  experimental 
observations. 

To  overcome  these  limitations,  we  have  developed  a  methodology  to  treat  recrystallization 
phenomena  using  electron  backscatter  diffraction  (EBSD)  measurements  and  two-dimensional  MC 
simulation  [5],  and  applied  these  techniques  to  model  the  recrystallization  of  cold-rolled 
commercial-purity  titanium.  Since  the  modeling  of  the  recrystallization  kinetics  of  work  hardened 
materials  has  important  industrial  applications  both  for  fabrication  and  for  the  prediction  of  in  service 
performance,  particular  attention  was  focused  on  the  effects  of  microstructural  parameters  on 
recrystallization  kinetics. 
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Experimental  procedures 

Commercial-purity  titanium  (cp-Ti  with  chemical  composition  of  600  O,  200  Fe,  15  H,  50  C,  100 
N,  balance  titanium  in  wppm)  was  chosen  as  the  program  material  for  this  investigation.  The 
as-received  material  was  cold  rolled  at  room  temperature  to  60%  total  thickness  reduction  and  then 
recrystallization  heat  treated  at  600°C  for  various  times. 

Microstructure  evolution  during  heat  treatment  was  established  using  EBSD  analysis  and 
micro-Vickers  hardness  measurements  were  conducted  to  determine  recrystallization  kinetics.  For 
EBSD  analysis,  samples  were  mechanically  polished  and  then  electro-polished  in  a  solution 
consisting  of  5  ml  perchloric  acid  and  95  ml  methanol  at  30V  and  -40°C.  For  EBSD  measurement,  a 
Hitachi  3400SE  field-emission-gun  scanning-electron  microscope  (FEG-SEM)  and  TSL-OIM™ 
software  were  used. 


Monte-Carlo  modeling  of  recrystallization 

The  details  of  MC  framework  and  recrystallization  model  employed  in  the  present  work  were  fully 
described  in  previous  work  [5],  and  thus  only  a  brief  description  is  provided  in  the  following. 
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Figure  1.  The  stored-energy  map  comprising  a  924x665 
triangular  lattice,  which  was  used  as  an  initial  input  for 
MC  modeling  of  recrystallization.  The  microstructure 
was  produced  by  mapping  the  EBSD  microstructure  of 
60%  cold-rolled  cp-Ti  onto  an  MC  grid.  The  step  size  of 
the  EBSD  scan  was  0.15  pm;  black  lines  represent 
high-angle  boundaries.  The  microstructure  contains 
2,167  grains  with  an  average  grain  size  of  30  pm.  Note 
that  the  fine  grains  show  higher  stored  energy,  compared 
to  the  coarse  grains. 


Construction  of  input  microstructure.  The  initial  microstructure  used  for  MC  modeling  was 
generated  by  mapping  the  EBSD  microstructure  of  60%  cold-rolled  cp-Ti  onto  the  MC  simulation 
grid  (Fig.  1).  Each  lattice  site  in  the  MC  grid  thus  generated  contains  information  on  the  orientation  in 
terms  of  Euler  angles  and  image  quality  (IQ)  value,  at  a  measurement  point  specified  by  the 
coordinate  pair  (x,  y)  from  an  area  measuring  120  pm  x  100  pm.  The  IQ  value  is  used  as  a  measure  of 
the  stored  energy  in  the  present  work,  and  the  stored  energy  of  each  lattice  site  is  calculated  as 
follows: 


Hi  =k(  1  - 


IQ, 

f-IQmax 


) 


(1) 


in  which  H;  and  IQ;  denote  the  relative  stored  energy  in  arbitrary  units  and  the  IQ  value,  respectively, 
of  the  MC  lattice  site  i,  and  IQmax  is  the  maximum  value  of  IQ.  The  factors  k  and  f  are  constants  which 
determine  the  upper  and  lower  bounds  for  the  stored-energy  distribution  and  were  set  as  70  and  1.3, 
respectively,  in  the  present  modeling. 


Recrystallization  model.  During  recrystallization,  the  microstructure  evolves  through  competition 
between  recovery,  nucleation  of  recrystallization,  and  grain  growth:  all  three  stages  of  the  process 
were  taken  into  account  in  the  present  modeling.  For  nucleation  kinetics,  site-saturated  nucleation 
was  assumed  so  that  all  nucleation  events  take  place  simultaneously  at  the  beginning  of 
recrystallization.  For  the  spatial  distribution  of  nuclei,  two  cases  were  considered:  random  nucleation 


2 


without  regard  to  location  in  the  microstructure  and  high-stored-energy  nucleation  (HSEN)  in  which 
nuclei  form  preferentially  at  high-stored-energy  regions  such  as  grain  boundaries,  triple  junctions, 
and  transition  bands.  The  movement  of  the  recrystallization  front  was  treated  with  the  strain-induced 
boundary  migration  (SIBM)  model  [6]  and  the  velocity  of  the  recrystallization  front,  Vy,  is  thus 
expressed  as  follows: 


Q 

■■  a-AHjj  -Mm  [  1  -  m-exp(  -  yjy  j3  7 


(2) 


where  a  is  proportionality  constant  and  m  is  a  parameter  which  controls  degree  of  anisotropic  grain 
boundary  mobility  [7];  AHy  and  0y  are  the  difference  in  stored  energy  and  misorientation  angle, 
respectively,  between  neighboring  sites  i  and  j;  Mm  is  the  maximum  mobility  corresponding  to  a 
random  high-angle  boundary.  The  velocity  of  the  recrystallization  front  in  Eq.  (2)  varies  from  0  to 
positive  values  less  than  unity,  depending  on  the  local  variation  of  the  stored  energy  and  the 
grain-boundary  mobility  throughout  the  microstructure.  When  vy  is  greater  than  a  random  probability 
obtained  from  a  random-number  generator  that  produces  numbers  between  0  and  1,  then  the 
recrystallized  grain  grows  at  the  expense  of  the  deformed  site  j;  as  a  consequence,  the  stored  energy  of 
site  j  is  reset  to  zero,  and  the  orientation  of  site  j  is  changed  to  that  of  site  i. 

The  recovery  process  during  recrystallization  was  simulated  such  that  when  a  chosen  deformed 
site  fails  to  undergo  recrystallization,  the  stored  energy  of  that  lattice  site  was  lowered  by  a  fixed 
factor,  b  (0<b<l),  without  changing  its  crystal  orientation.  Thus,  the  stored  energy  of  a  given  lattice 
site  which  underwent  recovery,  Hnew,  was  calculated  as  follows: 

Hnew  =  (l-b)Hold  (3) 

where  H0id  is  the  stored  energy  of  a  lattice  site  just  before  recovery  takes  place.  For  grain  growth, 
curvature-driven  grain  growth  was  assumed  [8]  as  is  in  a  conventional  MC  modeling. 


Results  and  discussion 

Recrystallization  kinetics  of  60%  cold-rolled  cp-Ti  at  600°C  was  determined  by  hardness 
measurements  and  EBSD  analyses.  The  kinetic  curve  exhibited  a  sigmoidal  shape  (Fig.  2(a))  which  is 
commonly  found  in  recrystallization  of  deformed  metals. 


Figure  2.  Recrystallization  kinetics  of  60%  cold-rolled  cp-Ti  during  heat  treatment  at  600°C:  (a)  plot 
of  fraction  recrystallized  against  logarithm  of  time,  (b)  the  JMAK  plot  of  (a),  and  (c)  the  SF  plot  of  (a). 
Note  that  experimentally-determined  RX  kinetics  deviates  from  linear  JMAK  kinetics  while 
exhibiting  linearity  when  plotted  based  on  the  SF  relation. 


The  experimental  data  were  re-plotted  in  Fig.  2(b),  based  on  Johnson-Mehl-Avrami-Kolmogorov 
(JMAK)  model  [9]  in  which  the  fraction  recrystallized,  Xrx,  depends  time,  t,  as  follows: 
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XRX=l-exp(-Btn) 


(4) 


in  which  B  is  a  parameter  that  depends  on  nucleation  and  growth  characteristics  and  nJMAK  is  the 
JMAK  exponent.  The  JMAK  kinetics  (Fig.  2(b))  exhibited  two  different  JMAK  exponents,  2.1  and 
0.7,  at  the  initial  and  later  stages  of  recrystallization,  respectively,  indicating  a  deviation  of 
recrystallization  behavior  of  the  present  material  from  the  JMAK  model.  On  the  other  hand,  the 
experimental  data  were  linear  when  plotted  based  on  the  empirical  Speich-Fisher  (SF)  parabolic 
relation  [10]: 


Xrx/(1-Xrx)  =  KtnsF  (5) 

in  which  K  is  a  constant  for  isothermal  reaction  and  nsF  is  the  SF  exponent.  Linearity  with  nsF  of  2.05 
was  found  in  the  SF  kinetics  (Fig.  2(c)),  thus  suggesting  that  the  recrystallization  kinetics  of  the 
current  material  can  be  well  explained  by  the  SF  relation  and  that  the  growth  rate,  G,  of  new  grain 
decreases  with  time,  i.e.,  G  is  proportional  to  1/t. 

In  order  to  investigate  the  effects  of  microstructural  factors  on  recrystallization  kinetics,  MC 
modeling  was  conducted  assuming  various  conditions  in  which  distribution  type  and  the  number  of 
nuclei,  stored  energy  distribution,  recovery  rate  and  grain  boundary  characteristics  were  varied;  the 
MC -predicted  recrystallization  kinetics  are  plotted  based  on  the  JMAK  relation  in  Fig.  3.. 


Figure  3.  MC-modeling  results  showing  the  effect  on  recrystallization  kinetics  of:  (a)  the  number  of 
nuclei,  (b)  average  stored  energy,  (c)  the  spatial  distribution  of  nuclei,  and  concurrent  recovery  under 
(d)  random  nucleation  and  isotropic  grain  boundary  motion,  (e)  HSEN  and  isotropic  grain  boundary 
motion  and  (f)  HSEN  and  anisotropic  grain  boundary  motion. 


It  was  found  from  the  modeling  that  increasing  the  number  of  nuclei  accelerates  the  overall 
recrystallization  kinetics  but  has  little  influence  on  the  slope  of  the  JMAK  plot  (Fig.  3(a)).  The 
increase  in  the  average  stored-energy  under  the  non-uniform-distribution  of  stored-energy  condition 
also  showed  similar  effects  to  the  number  of  nuclei  on  recrystallization  kinetics  (Fig.  3(b)).  A 
variation  of  the  distribution  pattern  of  nuclei,  on  the  other  hand,  was  found  to  have  a  considerable 
effect  on  the  JMAK  kinetics:  random  distribution  of  nuclei  resulted  in  a  linear  JMAK  kinetics  while 
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uniform  and  heterogeneous  distributions  of  nuclei  led  to  positive  and  negative  deviations, 
respectively,  from  the  theoretical  JMAK  kinetics  (Fig.  3(c)).  Concurrent  recovery  with 
recrystallization  also  affected  the  JMAK  kinetics  such  that  an  increase  in  recovery  rate  (i.e.,  increase 
in  the  parameter  ‘b’)  led  to  a  more  pronounced  deviation  from  a  linear  JMAK  kinetics  (Fig.  3(d)).  It 
was  also  revealed  that  a  combination  of  recovery  with  microstructural  inhomogeneities,  such  as 
heterogeneous  distribution  of  nuclei,  non-uniform  distribution  of  stored-energy,  and  anisotropic  grain 
boundary  mobility,  led  to  a  significant  deviation  from  the  theoretical  JMAK  kinetics  (Figs.  3(e)  and 
(f)). 


Figure  4.  Recrystallization  kinetics  predicted  by 
MC  modeling  under  two  extreme  conditions: 
one  satisfying  assumptions  in  the  JMAK  model 
and  the  other  including  the  effects  of 
heterogeneous  distribution  of  nuclei, 
non-uniform  stored-energy,  recovery  and 
anisotropic  energy  and  mobility  of  grain 
boundaries.  The  predicted  recrystallization 
kinetics  is  plotted  based  on  the  SF  relation  (Eq. 

(5))- 


Shown  in  Fig.  4  are  the  SF  plots  of  recrystallization  kinetics  which  were  predicted  by  MC 
modeling  under  two  extreme  conditions:  one  for  the  assumptions  in  the  classical  JMAK  model  and 
the  other  including  various  microstructural  effects  of  heterogeneous  distribution  nuclei,  non-uniform 
stored-energy,  recovery  and  anisotropic  energy  and  mobility  of  grain  boundary.  Recrystallization 
kinetics  under  the  former  assumptions  exhibited  a  positive  deviation  from  linearity,  reflecting  a 
constant  growth  rate  of  new  grains  in  the  JMAK  model.  On  the  other  hand,  recrystallization  kinetics 
predicted  using  the  latter  assumptions  showed  near-linearity,  which  is  in  good  agreement  with 
experimental  observation  (Fig.  2(c)). 

The  effect  of  such  microstructural  parameters  on  the  grain-size  distribution  was  also  predicted  by 
MC  modeling,  as  shown  in  Fig.  5.  The  experimental  results  showed  near  log-normal  distributions, 
which  were  unaffected  by  the  extent  of  recrystallization  (Fig.  5(a)).  When  effects  from  either 
non-uniform  stored-energy,  recovery  or  heterogeneous  distribution  of  nuclei  were  included  in  the 
JMAK  assumptions  (Fig.  5(b)),  MC  modeling  showed  that  the  former  two  effects  have  little  influence 
on  the  distribution  of  grain  size,  leading  to  a  narrow  distribution  skewed  toward  smaller  grain  sizes. 
On  the  other  hand,  a  broader  and  less  skewed  distribution  was  predicted  when  heterogeneous 
distribution  of  nuclei  was  included  in  the  JMAK  assumptions,  suggesting  a  critical  role  of  the 
distribution  pattern  of  nuclei  in  developing  a  log-normal  type  distribution.  When  all  of  the 
aforementioned  microstructural  effects  were  included,  MC  modeling  predicted  a  near  log-normal 
distribution  and  an  insensitivity  of  the  distribution  pattern  to  Xrx  (Fig.  5(d)),  both  of  which  are  in 
good  agreement  with  the  experimental  observations. 
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Figure  5.  Size  distributions  of  new  grains:  (a)  experimental  results,  and  the  distributions  predicted  by 
MC  modeling  under  conditions  (b)  that  include  an  effect  from  either  of  non-uniform  stored-energy, 
recovery  or  heterogeneous  distribution  of  nuclei  in  the  JMAK  assumptions  and  (c)  that  include  all  of 
the  effects  in  (b). 


Conclusions 

Both  experimental  observations  and  MC  modeling  of  recrystallization  of  cold-rolled  cp-Ti  reveal 
that  the  assumptions  in  the  classical  JMAK  model  result  in  a  failure  to  accurately  predict 
recrystallization  kinetics  while  the  SF  relation  provides  a  better  description.  MC  modeling  suggests 
that  the  inclusion  of  a  heterogeneous  distribution  of  nuclei  is  an  essential  prerequisite  for  more 
realistic  predictions  of  recrystallization  kinetics  and  grain-size  distributions. 
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